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ABSTRACT 


An experiment was performed to measure the scattering 
of underwater sound from a porous solid for the first time. 
Two porous solid spheres composed of heat-epoxied glass 
beads of 100 and 500 micron mean bead diameter were used. 
The permeability, porosity, and shear modulus of each sample 
were estimated from measurements made on cylindrical samples 
which had been manufactured at the same time and of the same 
glass beads as the spheres. These material properties were 
used as input to a theoretical model for the sound scattered 
from a poro-elastic sphere imbedded in a poro-elastic host 
developed by Kargl and Lim. The experimental data were 
compared to the theoretical calculations. 

Theoretical calculations with 0%, 3%, and 10% skeletal 
frame damping were compared to experimental data. Very good 
agreement between measured and predicted scattering was 
obtained for each sample over certain frequency ranges, 
taking 10% frame damping in the calculations. For other 
frequency ranges the agreement was less than good. No 
systematic trend in the agreement could be discerned with 


regard to porous grain size or sound frequency. 
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I. INTRODUCTION 


A. BACKGROUND 

There has been a resurgence of interest in sound 
propagation in a fluid-saturated porous medium because of 
the importance of ocean bottom propagation in a shallow 
water environment. While sound propagation in a porous 
medium is not a new subject, the scattering of elastic waves 
propagating in a porous medium has only recently been 
considered. An important potential application of this 
Scattering is the acoustic detection of buried mines. 

A general theory for the propagation of elastic waves in 
a porous, fluid-filled medium was derived by Maurice Biot 
(Biot, 1956a,b, 1962a,b). Biot’s theory has recently been 
applied to the scattering of elastic waves in a porous 
medium by a compact porous object (Kargl and Lim, 1993). To 
date, no experiments have been performed to validate this 


hallmark theoretical application of Biot’s theory. 


B. OBJECTIVES 

The objectives of this research are to measure and 
compare the scattering of underwater sound waves by a fluid- 
Saturated poro-elastic sphere to theoretical values computed 


following Kargl and Lim (Kargl and Lim, 1993). 


It is noteworthy that this experiment, the scattering of 
underwater sound from porous solid spheres of glass beads, 


is the first of its kind. 


II. HIGHLIGHTS OF SOUND PROPAGATION AND SCATTERING IN 


POROUS MEDIA 


In a series of papers, Maurice Biot (Biot 1956a,b, 
1962a,b) developed a theory for the propagation of elastic 
waves in a fluid-saturated porous solid.  Biot derived 
equations that describe the interaction between a fluid- 
saturated porous solid and the fluid in the pores due to the 
solid's microscopic pore geometry. 

The essence of Biot theory is that among the acoustic 
properties in a fluid-saturated porous solid there exist 
three propagating body waves: two longitudinal, so called 
fast and slow compressional waves, and a shear wave. The 
existence of the slow wave was first confirmed by Plona 
(Plona, 1980, pp. 259-261). Other measurements have also 
confirmed the existence of the slow wave (Mayes et. al., 
1986, p.249). The slow wave is characterized by 1) the 
sense of the motion of the fluid and solid are opposite, 
i.e. "out of phase", 2) it travels with a speed that is 
Slower than the bulk speed of the fluid, and 3) it is 
typically very lossy. The fast wave is characterized by 1) 
the sense of the fluid and solid motion are the same, i.e. 
"in phase", 2) the speed of this wave is between that of the 


bulk speed of the fluid and the solid, and 3) the 


attenuation is small. When we talk about sound propagation 
in the ocean bottom we are talking about the fast wave. The 
shear wave can be described as a shear wave that propagates 
in the porous solid frame with the fluid dragged along by 
viscous and inertial interaction. 

The speed and attenuation of these waves are determined 
by eleven fundamental material parameters: three material 
properties of the bulk fluid: 1) mass density, p, 2) 9p 
modulus, K, and 3) shear viscosity, 7, five material 
properties of the solid scattering object: 1) porosity, P, 
2) mass density, p,, 3) bulk modulus of the solid grains, K, 
4) bulk modulus of the sketal frame, K,, and 5) shear 
modulus of the skeletal frame, N, and three parameters that 
depend upon the microscopic pore geometry: 1) tortuosity, o, 
2) permeability, k, and 3) structural factor, ô, in ENS 


universal frequency-dependent correction function. 


ІІІ. SAMPLE PREPARATION AND MATERIAL PROPERTY DETERMINATION 


A. DESCRIPTION OF SAMPLES 

The porous spheres used as scatters in the experiment 
are composed of bonded borosilicate glass beads of nominally 
the same diameter. The unique manner in which the spheres 
were fabricated is worth discussion. 

Many different types of experiments have used 
consolidated glass beads to simulate sediment (Johnson and 
Plona, 1982, p. 557). The usua. process by which glass 
beads are consolidated is heat fusion. This method was 
initially attempted to fabricate the glass bead spheres for 
the present experiments but proved to be inadeguate. 

The glass beads were poured into a spherical mold and 
heated to a temperature of approximately 704 C. Тһе 
softening temperature of borosilicate is 821 C (Moses, p. 
869). The samples were "cooked" for 12 hours. After a 
cooling period the samples were inspected. It was found 
that the glass beads had settled, resulting in a sphere-like 
sample with a smooth spherical lower hemisphere but with a 
squashed upper hemisphere. 

Different temperature and cook-time combinations were 
investigated in an attempted to reduce the settling. When 


the cook time was decreased to 8 versus 12 hours, or the 


cook temperature reduced a significant amount to 593 C 
versus 704 C, the degree of settling was reduced. However, 
the glass beads scarcely adhered to one another, producing 
very fragile spheres with beads that would crumble off at 
the slightest touch. When the cook time or cook temperature 
was increased, the amount of settling increased. There was 
no cook time and temperature combination that produced a 
Satisfactory sphere. 

It was decided to have the samples commercially 
manufactured by a company experienced in fabricating bonded 
glass beads (General Polymeric, Reading PA). To make the 
samples, the glass beads are coated with a fine powder 
epoxy, then heat fused at a temperature sufficient to 
activate the epoxy, 150 C. Settling of the glass beads was 
no problem because the curing temperature was much lower 
than the softening temperature of the glass beads. 

Borosilicate, general purpose, clear glass beads were 
used to manufacture two spheres. The same mold was used for 
both. One sphere was composed of glass beads with a mean 
bead diameter of 100 microns while the other of glass beads 
with a mean bead diameter of 500 microns. Both spheres had 
a major and minor axis diameter of 0.066 meters. 

Due to the difficulty in measuring the material 
properties of spherical samples, cylindrical samples of each 


bead size were prepared at the same time as the spheres. 


These samples were 0.151 meters long and 0.052 meters in 


diameter. 


B. DETERMINATION OF POROSITY 
The ratio of the open pore volume, the volume not 
occupied by glass, to the total volume of the sample (sphere 


or cylinder) is its porosity, P. 


open pore volume (3.1) 


P = 
total volume 


The porosity of the cylinders was determined as the 
ratio of the weight of the cylindrical porous sample to the 
calculated weight of the sample if were solid borosilicate 


glass. The porosity of each is listed in Table 3.1. 


TABLE 3.1 POROSITIES OF THE TWO CYLINDERS AND THE VALUES 


USED IN THEIR CALCULATIONS 
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Density of Glass 2231 kg/m 2231 kg/m 















Beads, 0. 


Calculated Mass, M, 0.738 kg 0.730 kg 
— 


С. DETERMINATION OF PERMEABILITY 







Permeability, k, is a measure of the ease with which a 
porous solid allows fluid to flow through it. It is 
proportional to the square of the nominal pore radius. 
Permeability is related to the DC flow resistance of a fluid 
saturated porous solid. DC flow resistance, Rpc, is the 
average gradient in fluid pressure, Vp, in the porous solid 
divided by the average fluid velocity, dU/dt in the 
direction of the gradient (Baker, 1986,p. 34), and is given 


in Equation 3,27 


EDENE (3.2) 


In terms of the permeability k, 


ee md (3.3) 


where P is the porosity and y is the fluid shear viscosity. 
Note that, unlike the DC flow resistance, the permeability 
is a property of the microscopic geometry of the porous 
solid only, independent of the filling fluid. 

The permeability for each sphere was experimentally 
determined by the constant head permeability test (Means and 
Parcher, 1986, pp. 407-417). In the constant head 
permeability test the volume of water that passes through 
the porous sample over a period of time, along with the 
sample's dimensions, is used to calculate a coefficient of 
permeability, k’, with units of meters per second (The 
symbol k’ is being used instead of k for Means and Parcher’s 
coefficient of permeability to avoid any confusion with 
permeability defined previously) 


Figure 3.1, adapted from Means and Parcher, Physical 


Properties of Soils, Charles E. Merrill Inc., Columbus, Ohio 
(1964), illustrates the constant head permeameter used to 


collect the data. 








To Water 
Reservoir 


Adjustable 
Mounting 


e Evaporating 
Dish 


Sink 


Figure 3.1. Constant Head Permeameter 


The cylindrical samples of length L and diameter d were 
securely wrapped in "plumber's tape" and connected to one 
end of a lucite tube using a rubber gasket. The entire 
apparatus was attached to a lab stand over a sink. A 
constant supply of tap water was applied to the top of the 


lucite tube. Once the desired water column height was 


10 


established, a constant head was maintain through the use of 
an overflow tube. After a sufficient period of time had 
elapsed, ensuring the sample was thoroughly saturated, the 
volume of water that passed through the sample over a period 
of time was measured. 

Equation 3.4 is used to compute the coefficient of 


permeability, xk” 
pri QL (3.4) 


where Q is quantity of water collected, L is sample length, 
t is time required to collect the water, A is area of the 
sample and h is the distance between the bottom of the 
sample and the height of the water in the lucite tube. 

То get a value of permeability in mks units (m’), the 
coefficient of permeability must be multiplied by the 


conversion factor 
aL 71.02 x 107 (3.5) 


where y is the shear viscosity of water (0.001 Pas), g 15 
the acceleration of gravity and p is the density of the 


water (998 kg/m) (Baker, 1986, P. 33). Table 3.2 contains 


Mur 


the experimentally determined values used to compute the 
coefficient of permeability, k', of each sample along with 


the permeability, k. 


TABLE 3.2. CONSTANT HEAD PERMEABILITY DATA 


AA 














D. DETERMINATION OF ELASTIC MODULI 

The last two parameters determined experimentally were 
the shear and bulk moduli. The resonant acoustic method was 
used to determine the elastic moduli of the cylindrical 


samples (Garrett, 1990). The resonant acoustic method is 


то 


based on the notion that a uniform cylinder, with a cross 
section of diameter d and of length L which is significantly 
greater than the diameter, will propagate three independent 
waves (longitudinal, torsional, and flexural) if their 
wavelengths A are much greater than the diameter d. These 
modes will exhibit resonance at frequencies that depend upon 
the length and the boundary conditions at the ends of the 
rods. 

In the classic resonant acoustic experiment a long 
cylindrical rod is selectively excited in its torsional, 
longitudinal, and flexural modes using a single pair of 
electrodynamic transducers consisting of coils glued to each 
end of the bar and placed in the field of a magnet. See 
Garrett for a through discussion of the experiment (Garrett, 
NOS). 

The previously introduced porous cylindrical rods were 
prepared for this measurement in a similar fashion. The 
transducer coils were 20 windings of magnet wire, 
approximately one meter in length, epoxied to each end of 
the rods. Particular attention was paid to limiting the 
amount of epoxy used so as to minimize the mass loading of 
the bar ends by the transducers. The total mass of each 
pair of coils and the epoxy used to adhere it to the rod was 
5.15 grams. Since the mass of the 100 micron sample was 


510.41 grams and the mass of the 500 micron sample was 
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495.23 grams, any mass loading caused by the coils and epoxy 
is considered negligible. 

The optimum sample aspect ratio, L:d, for the acoustic 
resonant experiment is on the order of 20:1 (Garrett, 1990, 
p. 212). Тһе two porous cylinder samples have aspects 
ratios of approximately 3:1. Ав a result, several apparatus 
modifications were necessary. 

The major impact of the relatively small aspect ratio 
was excessive electrical crosstalk between the two end coils 
due to magnetic flux coupling. Two modifications were made 
to reduce the crosstalk as much as possible. The first 
modification was to the alignment of the cylindrical rod's 
two transducer coils. Normally the coils are oriented with 
their planes parallel to one another. To reduce the coupled 
magnetic flux, the planes of the coils were instead oriented 
90 degrees with respect to each other. The second 
modification was to wrap a moveable "humbucking" coil around 
the middle of the cylinder. The position of and drive level 
to this coil were adjusted to null out as much of the 
crosstalk as possible. 

Each cylindrical rod was placed with its axis at the 
center of the pole-piece faces, which were aligned along the 
bar axis. Figure 3.2, from Garrett, "Resonant Acoustic 
Determination of Elastic Moduli," J. Acoust. Soc. Am. 88, 


213 (1990), shows the magnet/rod configuration at each end. 
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Figure 3.3 illustrates the positioning of the coils with 


respect to the magnet pole pieces. 





Figure 3.2. Position of the Transducer Coils Relative 
to the Magnet's Pole Pieces for the Torsional Modes. 





Figure 3.3. Modified Rod and Magnet Configuration for 
the Acoustic Resonant Experiment. 
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The acoustic resonant experiment can be used to 
determine all three resonant modes, torsional, flexural and 
longitudinal, of a long thin sample rod. Although attempts 
were made to observe all three modes, due to the excessive 
electrical crosstalk, only the torsional mode was 
observable. For both the 100 micron and 500 micron samples 
only the first and second — modes were observable. 
For these it was possible to measure only their resonance 
frequency, not their quality factor. 

The first and second harmonic modes were observed in 


both the 100 micron and 500 micron samples. The shear 


modulus N was computed using equation 3.6: 


N = ApL? (f, /n)? (3.6) 


where L is the length, f, is the frequency of the nth mode, 
n is the mode number, and p is the mass density of the rods. 


p is computed using equation 3.7: 


procu D (3.7) 


Where p is the bulk mass density of the glass beads. The 
shear modulus for each sample is listed in Table 3.3. 


Since the longitudinal modes were not observable, the 
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bulk modulus had to be estimated using the shear modulus and 
Poisson's ratio. Equation 3.8 shows the relationship 
between the Poisson's ratio, v, and the shear, N, and bulk, 


K, modulus of a material. 


у = .3iK -2N (3.8) 
23K FN) 


The exact value of Poisson's ratio for the sample rods 
is not known.  Poisson's ration for the borosilicate glass, 
of which the beads are made, is 0.20 (Moses, p. 869). There 
was some concern as to whether Poisson's ratio for a porous 
rod composed of borosilicate beads should be the same as 
bulk borosilicate. It was noted that Berryman's effective 
medium calculation of the elastic moduli for the porous 
solid samples, which is based upon input elastic moduli for 
borosilicate, and which is used by Kargl and Lim, yields a 
Poisson's ratio of 0.20 for the porous sample. Accordingly, 
a Poisson's ratio of 0.20 was assumed in calculating the 
bulk modulus for the porous solid samples from the measured 
shear modulus. The bulk and shear moduli for each sample 


are listed in Table 3.2. 
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TABLE 3.2.  POROUS SAMPLE SHEAR AND BULK MODULI 


14° 07x10 Pa 13.62x10° Pa 























Shear Modulus, N 


Average Shear 


Bulk Modulus, K 













assuming v=0.20 
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IV. BACKSCATTER EXPERIMENT 


A. EXPERIMENT OBJECTIVE 
The backscattering experiment involved measuring the 
amplitude of the backscattered acoustic pressure from a 


porous solid sphere in water. 


B. EXPERIMENT COMPONENTS 

The backscatter experiment consisted of four major 
components: 1)a water-filled tank, 2)scattering objects 
(porous spheres), 3)an acoustic projector, and 4)an acoustic 
receiver (hydrophone). The experiment was conducted ina 
water-filled laboratory tank that was 7.32 meters in length, 
1.64 meters in width and 2.02 meters in depth. All tank 
boundaries, except the top surface, were covered with sound 
absorbing tiles. The scattering objects were the previously 
discussed porous spheres 6.64 centimeters in diameter. 

The acoustic projector was a type F33 general purpose 
directional transducer composed of two concentric arrays of 
piezoelectric vibrators (USRD Transducer Catalog 1991). 
Figure 4.1 is an illustration, from the USRD Transducer 
Catalog, April 1991, of the type F33 transducer. The inner 
array consists of 12 PZT disks that form a 3.8 by 5.0 


centimeter rectangle. The frequency range of the inner 
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array is 15-150 kHz. The outer array consists of 64 PZT 
disks that form a 20.8 by 21.7 centimeter rectangle. The 


frequency range of the outer array is 1-150 kHz. 





Units are cm 


Figure 4.1 Type F33 General Purpose Directional 
Transducer. 


The acoustic receiver was a Bruel and Kjaer type 8103 
miniature hydrophone. The type 8103 is illustrated, from 
its calibration chart, in Figure 4.2. It has an operating 


range of 0.1-180 kHz, is 2.5 centimeters long, has a radius 


of 0.475 centimeters. 
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Units are mm 


Figure 4.2 Bruel & Kjaer Type 8103 Miniature 
Hydrophone. 


C. LABORATORY TANK SETUP 

In order to compare the backscatter data to Kargl and 
Lim's theory, it was necessary to ensure that far-field 
conditions were met for the scatterer and the hydrophone 
(Kinsler et. al., p. 188). Equation 4.1 was used to compute 


the far-field distance, г, 


ee (4.1) 
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where d is the source diameter and A is the source 
wavelength. The values computed for radiation from the 
inner and outer arrays of the F33, and from the spherical 
scatters, are listed in Table 4.1 for the worst case 


frequency of 150 kHz. 
TABLE 4.1 FAR-FIELD DISTANCE 
Component Spheres Inner Outer 
Array Array 


Upon examination of Table 4.1 it is clear that the far- 













field distance for the F33 outer array exceeds the 
dimensions of the tank. For this reason the inner array was 
used as the projector. The components, projector, sphere 
and hydrophone were arranged in a configuration that 
Satisfied all far-field criteria. 

Figure 4.3 is an illustration of the tank and components 
configuration. The projector was position at one end of the 
tank, 1.01 meters below the surface, 0.82 meters from the 


sides, facing down the center of the tank. Prior to being 
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placed in the tank, the sphere being investigated was 
secured in a hairnet, placed in a beaker full of water from 
the tank, and vacuum degassed. This was necessary to ensure 
the sphere was as completely water-saturated as possible. 
The sphere was then lowered into the tank, still in the 
beaker full of water. It was then suspended from a nylon 
string attached to the hairnet 1.01 meters below the 
surface, on axis with the projector, 1.14 meters from the 
projector. The miniature hydrophone was the only component 
that did not remain at one location throughout the 
experiment. 

The hydrophone was positioned by suspending it from a 
"hula-hoop", which was secured, concentric with the sphere, 
to the grating that covered the tank. The radius of the 
"hula-hoop" was 0.44 meters and was marked at 10 degree 
increments to aid in the positioning of the hydrophone. A 
weight was attached to the lower end of the hydrophone cable 
by a length of nylon string, with a length such that the 
weight was positioned centimeters above the tank bottom. 

The miniature hydrophone was initially position on axis 
between the projector and the sphere a distance of 0.44 
meters from the sphere. This distance was chosen to ensure 
that the hydrophone was in the sphere's as well as the 
projector's far field. During the experiment the hydrophone 


was moved around the sphere in 10 degree increments. 
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Side View 


Hydrophone 
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Figure 4.3 Top and Side View of Laboratory 
Tank Setup. 


24 


D. ELECTRONICS 

Figure 4.4 is a block diagram of the electronics used in 
the experiment. The transmitted signal was generated by a 
Hewlett-Packard Model 3314A Function Generator and amplified 
by a Hewlett-Packard Model 467A Power Amplifier, then 
applied to the type F33 transducer. The maximum output 
voltage the amplifier is 20 volts peak-to-peak. The output 
of the signal generator was monitored with a Nicolet Model 
Pro 30 Digital Oscilloscope. The output from the Bruel and 
Kjaer type 8103 miniature hydrophone was filtered by a 
Ithaco 1201 Low Noise Preamplifier and input to the Nicolet 
Digital Oscilloscope for analysis. An Iwatsu SS5710 
Oscilloscope was used to monitor the waveforms from the 
power amplifier (transmitted signal) and the output from the 


low noise preamplifier (received signal). 
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HP 3314A HP 4674 
Function Generator Power Amplifier 
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Preamplifier 


Nicolet Digital! 
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Figure 4.4 Diagram of Electrical Components. 


E. EXPERIMENTAL PROCEDURE 

Based upon the frequency range of the inner array of the 
type F33 projector (USRD Catalog 1991), a frequency range of 
30 kHz to 150 Khz was used for the backscatter measurements. 


Backscatter measurements were taken from 150 kHz down to 
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30 kHz in 10kHz intervals for the 100 micron sphere and in 
30 kHz intervals for the 500 micron sphere. The relatively 
close proximity of the sphere, projector, hydrophone, and 
tank boundaries to each other determined the number of 
cycles (8) that was chosen for the transmitted pulse. 
Separating the signal of interest from the many boundary 
reflected signals was difficult at, and impossible beyond, 
angles of + 090 degrees from the 180 degree backscatter 
direction. 

With the hydrophone at its initial position, on axis 
between the projector and the sphere, the function generator 
was first set to 150 kHz. The typical voltage applied to 
the projector was 20 volts peak-to-peak. The filtered 
scattered signal from the hydrophone was averaged 512 times 
by the Nicolet digital oscilloscope. Figures 4.5 and 4.6 
are typical averaged screen displays from the Nicolet 
digital osccilloscope. The positive and negative peak 
values from the averaged signal were recorded, averaged, 
then added to produce an average peak-to-peak voltage of the 
Scattered return for that particular frequency/bearing 
combination. Particular care was taken to ensure that the 
ring-up and ring-down of the signal was not measured as part 
of the scattered return signal. The frequency of the signal 
was then decreased either 10 kHz or 30 kHz, depending upon 


which sphere was being investigated, to the next lower 
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De И 
Signal Screen Dump 





Figure 4.6 Typical Averaged Backscatter Signal Screen Dump 
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frequency and the scattered return measured. This procedure 
was continued down to a frequency of 30 kHz. The hydrophone 
was then moved in 10 degrees increments to the right and the 
procedure repeated until the data were no longer readable 
due to boundary interference. The left side was 
investigated in the same manner. After all the data were 
collected for the sphere it was lowered to the tank bottom 
and the hydrophone was positioned in the sphere's normal 
position. Next, the peak-to-peak values of the pulse 
incident at the hydrophone were collected over the same 
frequency range as for the scattering measurements. These 
peak to peak voltages were used in the normalization of the 
Scattered signals. The peak to peak scattered and incident 
values for the 100 micron sphere are listed in Appendix A, 
while the values for the 500 micron sphere are listed in 


Appendix B. 
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V. DATA ANALYSIS AND COMPARISON 


A.  EXPERIMENTAL RESULTS 
For comparison purposes, the peak-to-peak backscattered 
hydrophone output voltage values were normalized by the 


following equation: 


V 
Normalized Scattering = “= — (5.1) 


where V, is the peak-to-peak voltage of the incident 


Signal, V, 


5 


са 19 the peak-to-peak voltage of the scattered 
Signal, and r is the distance from the scatterer to the 
receiver, 0.44 m. Spherical spreading of the scattered wave 
was assumed. The values of the normalized scattering for 
all data taken are given in Appendix A and B, respectively, 


for the 100 micron and 500 micron samples. 


B. SOURCES OF THEORETICAL VALUES 

Theoretical values of scattered pressure were computed 
for a rigid and massive sphere following Morse (Morse, p. 
355) and for a porous sphere using a Fortran program 
provided by Kargl (Kargl, 1993) based upon Kargl and Lim’s 


application of Biot’s theory (Kargl and Lim, 1993). 
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Kargl's program calculates the scattering of fast and 
slow compressional waves from a fluid-saturated poroelastic 
sphere embedded in a fluid-saturated poroelastic host. Due 
to the versatility of the program, 29 inputs are required. 
In order to compare results from Kargl's program to the 
laboratory data for a porous sphere in water, the input 
parameters for the host porous medium were chosen to 
simulate water as closely as possible. The input values 
used are listed in Appendix C. 

As a confidence test of Kargl's program, input 
parameters were chosen to simulate a rigid and massive solid 
sphere target in water for comparison to known results. A 
rigid and massive sphere was approximated by a sphere with a 
density 100 times that of water, with a bulk modulus 10,000 
times that of water, and a Poisson's ratio of 1/3. The 
results of Kargl's program were compared to the scattering 
of a infinitely rigid and massive sphere in water following 
Morse (Morse, p. 355) for the frequencies of interest, 30, 
60, 90, 120 and 150 kHz. The comparisons are illustrated in 
Figures 5.1-5. The solid curves correspond to the 
theoretical values generated by Kargl's program while the 
dotted curves correspond to theoretical values generated 
following Morse. The agreement is satisfactory and is taken 


as validation of Kargl's program. 
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Figure 5.1. Normalized Scattering, 30 kHz Solid Sphere: 
Morse(o) Kargl(-) 





Figure 5.2. Normalized Scattering, 60 kHz Solid Sphere: 
Morse(o) Kargl(-) 
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Figure 5.5. Normalized Scattering, 150 kHz Solid Sphere: 
Morse(dotted) and Kargl(solid) 


C. COMPARISON OF EXPERIMENTAL RESULTS AND CALCULATIONS 
Figures 5.6 through 5.15 show a comparison of the 
measured scattering with theoretical calculations for the 
100 and the 500 micron spheres at five frequencies, 30, 60, 
90, 120, and 150 kHz. The solid curve in each figure was 
generated by Kargl’s program using the input data listed in 
Appendix C. The dashed curve is the scattering from a rigid 
and massive sphere following Morse. The open circles 
connected by dot-dash are data, where the dot-dash between 
the data points is simply a guide to the eye. It is obvious 
from the figures that the experimental and theoretical 
values do not agree for the input parameters used in the 


calculations. 
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gure 5.6. Normalized Scatterng, 30 kHz 100 Micron: 
Kargl(solid),Morse(dashed), & Exp. Data(dots superimposed) 





Figure 5.7. Normalized Scattering, 60 kHz 100 Micron: 
Kargl(solid), Morse(dashed), & Exp. Data(dots superimposed) 
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igure 5.8. Normalized Scattering,90 kHz 100 Micron: 
Kargl(solid), Morse(dashed), & Exp. Data(dots superimposed) 





ized Scattering, 120 kHz 100 Micron: 


Figure 5.9. Normal 
Kargl(solid),Morse(dashed), & Exp. Data(dots superimposed) 
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Figure 5.10. Normalized Scattering, 150 kHz 100 Micron: 
Kargl(solid), Morse(dashed), & Exp. Data(dots superimposed) 





Figure 5.11 Normalized Scattering, 30 kHz 500 Micron: 
Kargl (solid), Morse (dashed) ,‚& Exp. Data (dotssuperimposed) 
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Figure 5.12 Normalized Scattering, 60 kHz 500 Micron: 
Kargl(solid), Morse(dashed), & Exp. Data(dots superimposed) 





Figure 5.13 Normalized Scattering, 90 kHz 500 Micron: 
Kargl (solid), Morse(dashed), & Exp. Data(dots superimposed) 
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Figure 5.14 Normalized Scattering, 120 kHz 500 Micron: 
Kargl(solid), Morse(dashed), & Exp. Data(dots superimposed) 





Figure 5.15 Normalized Scattering, 150 kHz 500 Micron: 
Kargl (solid), Morse(dashed), & Exp. Data(dots superimposed) 
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The most notable difference between the experimental апа 
theoretical results are in (1) the overall amplitude of the 
scattering and (2) the structure in the angular dependence 
of the scattering. Generally, the observed scattering is 
much less than predicted. Regarding the angular dependence, 
the observed results resemble much more the scattering from 
a rigid and massive sphere with reduced amplitude, than the 
predicted results. It is not likely that the lack of 
Structure in the observed scattering is a consequence of 
sampling the angular function every 10 degrees. In numerous 
cases there are several experimental samples within the 
angular ranges of a lobe in the predicted results. 
Furthermore, the 180 degree backscatter was measured at all 
frequencies, and no hint of strong lobing was experimentally 


observed in this direction. 


D. SPECULATION FOR DISAGREEMENT BETWEEN THEORETICAL 

AND EXPERIMENTAL VALUES 

The amplitude and structure of the calculated 

predictions for the scattering are suggestive of internal 
resonance in the scatter. It was speculated that the 
presence of these features may be due to the neglect of 
Material damping of the porous solid frame in the calculated 
scattering. To test this speculation two additional 


calculations of scattering using Kargl’s program were 
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- 


performed. In the first case a damping value of 3$ was 
used. This value was chosen based on Brown's measurement of 
the elastic modulus of a common epoxy (Brown, 1991, p. 41). 
The second test was performed using a damping value of 10$. 
This represent's a very large damping value. The values of 
the elastic moduli use in the calculations corresponding to 
these damping values are listed in Appendix C. 

Figures 5.16 through 5.25 are comparisons between the 
theoretical values produced by Kargl's program with 3$ 
damping (dashed curve), 10$ damping (dotted curve), and 


experiment data (solid curve with dots superimposed). 





Figure 5.16 Norm. Scat., 30 kHz 100 Micron: 
3% Damping (dashed), 10% Damping (dotted), 
& Experiment Data(dots superimposed) 
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Figure 5.17 Norm. Scat., 60 kHz 100 Micron: 
3% Damping (dashed), 10% Damping (dotted), 
& Experiment Data(dots superimposed) 





Figure 5.18 Norm. Scat., 90 kHz 100 Micron: 
3% Damping(dashed), 10% Damping (dotted), 
& Experiment Data(dots superimposed) 
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Figure 5.19 Norm. Scat., 120 kHz 100 Micron: 
3% Damping (dashed), 10% Damping (dotted), 
& Experiment Data (dots superimposed) 





Figure 5.20 Norm. Scat., 150 kHz 100 Micron: 
3% Damping(dashed), 10% Damping(dotted), 
& Experiment Data(dots superimposed) 
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Figure 5.21 Norm. Scat., 30 kHz 500 Micron: 
3% Damping (dashed), 10% Damping (dotted), 
& Experiment Data (dots superimposed) 





Figure 5.22 Norm. Scat., 60 kHz 500 Micron: 
3% Damping (dashed), 10% Damping (dotted), 
& Experiment Data (dots superimposed) 
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Figure 5.23 Norm. Scat., 90 kHz 500 Micron: 
3% Damping (dashed), 10% Damping (dotted), 
& Experiment Data(dots superimposed) 





Figure 5.24 Norm. Scat., 120 kHz 500 Micron: 
3% Damping(dshed), 10% Damping(dotted), 
& Experiment Data(dots superimposed) 
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Figure 5.25 Norm. Scat., 150 kHz 500 Micron: 
3% Damping (dashed), 10% Damping (dotted), 
& Experiment Data(dots superimposed) 


From the figures it is apparent that including material 
damping in the skeletal frame substantially improves the 
agreement between measured and calculated scattering. Very 
close agreement is found between the experimental data and 
Kargl's curves for 10% damping for the 500 micron sample at 
frequencies of 90 kHz and above (Figures 5.23 through 5.25). 
Except for the pronounced lobe in the 180 degree backscatter 
direction predicted by the calculations, the agreement 
between the data and the 10% damping curves is fair for the 
100 micron sample between 30 and 90 kHz (Figures 5.16 


through 5.18). In the remaining cases the agreement is less 
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than good, with no recognizable trend with frequency or 
grain size (Figures 5.19 through 5.22). 

As mentioned previously, a damping value of 10$ is very 
large. It is a legitimate question as to whether this is 
realistic, even though the agreement between the data and 
10$ damping curves was better than for the 3$ damping 
curves. As was mentioned previously, if greater aspect 
ratio (20:1) porous cylinder samples were available, the 


Skeletal damping could be directly measured. 
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VI. SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS 


A. SUMMARY AND CONCLUSIONS 


An experiment was performed to measure the scattering of 
underwater sound from a porous solid for the first time. 

Two porous solid spheres composed of heat-epoxied glass 
beads of 100 and 500 micron mean bead diameter were used. 
The permeability, porosity, and shear modulus of each sample 
were estimated from measurements made on cylindrical samples 
which had been manufactured at the same time and of the same 
glass beads as the spheres. These material properties were 
used as input to a theoretical model for the sound scattered 
from a poro-elastic sphere imbedded in a poro-elastic host 
developed by Kargl and Lim. The experimental data were 
compared to the theoretical calculations. 

Theoretical calculations with 0%, 3%, and 10% skeletal 
frame damping were compared to experimental data. Very good 
agreement between measured and predicted scattering was 
obtained for each sample over certain frequency ranges, 
taking 10% frame damping in the calculations. For other 
frequency ranges the agreement was less than good. No 
systematic trend in the agreement could be discerned with 


regard to porous grain size or sound frequency. 
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В. RECOMMENDATIONS 

Regarding the scattering, the greatest hinderance to 
obtaining better quality data was the small size of the 
tanks available at NPS. It is recommended that future 
measurements be conducted in a large tank. This would allow 
Scattering measurements to be made down to lower 
frequencies, where the viscous drag on the fluid in the 
pores is more significant, and would also allow more of the 
angular dependence to be measured. 

Regarding the calculation of the scattering, it is 
believed that all the material properties of the porous 
Spheres are fairly well known except the sketal frame 
damping. Since it has a great influence on the scattering, 
it would be desirable to measure this directly. If greater 
aspect ratio (20:1) porous cylinder samples were available, 
both the torsional and longitudinal moduli and the damping 
coefficients could be measured. 

Lastly, no attempt was made to measure the tortuosity of 
the present samples; a typical value was assumed based upon 
the porosity. Although it would not be expected to differ 


significantly, the tortuosity could be measured. 
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Appendix A: Experiment Data 100 Micron Sphere 


Scattered Incident Ratio of Radius, Normalized 
Voltage Peak Voltage Peak | Scattered & Voltage 


— | 





90 Degrees 
Frequency 






















| 1 0.0135 
30K Hz | 7.1300 | 0.1545 0.0461 0.4450 | 0.0205 






100 Degrees | SCO 






150 K Hz 47.2183 6.4486 0.0073 0.4450 0.0033 
140 K Hz 48.3593 5.6726 | 0.0085 0.4450 0.0038 







130 K Hz 48.7275 51975 0.0091 0.4450 0.0040 
120 K Hz 46.1983 5.2122 0.0089 0.4450 0.0039 


















110 K Hz 45.7013 4.8513 | 0.0094 0.4450 0.0042 
100 K Hz 41.0225 3.6166 | 0.0113 0.4450 0.0050 
90K Hz | 19.5000 | 2.3823 | 0.0082 0.4450 0.0036 









80K Hz | 17.7602 1.5898 | 0.0112 0.4450 0.0050 
70K Hz | 17.2715 i ç 0.4450 0.0065 





17.6110 : 0.0091 
50 K Hz | 11.0512 0.5325 | 0.0208 0.4450 0.0092 








40 K Hz 8.5500 0.3127 | 0.0273 0.4450 0.0122 


30K Hz | 9.5092 0.1545 | 0.0615 0.4450 0.0274 
| 


| 
o Degrees | T | | ши 
Гвокнеј 174180 | 1588 | 00110 10.4450] 00049 | 
30 K Hz | 1.4373 0.1545 
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Appendix A: Experimental Data 100 Micron Sphere 


120 Degrees | | 


































Frequency Scattered Incident Ratio of Radiu d Normalized 
Voltage Peak Voltage Peak Scattered & Voltage 
120K Hz | 415385 
33.2975 | 2.3823 0.0140 
80 K Hz 17.2008 1.5898 0.0108 
70 K Hz 16.5760 | 1.2265 | 0.0135 
50K Hz 4.5138 | 0.5325 | 0.0085 
3.9103 | 0.3127 | 0.0125 0.4450 
30 K Hz 1.0550 | 0.1545 | 0.0068 | 0.4450 
ЕЕ | | | | 
130 Degrees | | | 
150 K Hz 415.0708 — — 6.4486 0.0070 | 0.4450 0.0031 
140 K Hz 39.5708 —— 5.6726 | 0.0070 0.4450 
394028 — — | 0.4450 0.0033 
437050 | 1193 | 0.0084 
52.1192 | 4.8513 | 0.0107 
136108 | К — 
90 K Hz 23.3100 2.3823 | 0.0098 
80 K Hz 226865 | 1.5898 | 0.0143 
70 K Hz 12.7390 /| 1.2265 | 0.0104 0.4450 0.0046 
60K Hz 11.4392 | 0.8601 | 0.0133 
50K Hz 5.5797 0.5325 | 0.0105 
40 K Hz o | 0.3127 | 0.0173 
30KĦz | 3225 01545 | 0.0211 





| 
22.4520 1.5898 | 0.0141 
"WES | 05325 0.0132 
40 K Hz 5.4743 0.3127 
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Appendix A: Experimental Data 100 Micron Sphere 


150 Degrees ар (ии 

Frequency Scattered Incident | Ratio of Normalized 
Voltage Peak Voltage Peak Scattered & Voltage 

| 60KHz | 102267 | 08601 | 0019 (0445 0.0053 | 

4.0188 0.1545 | 0.0260 0.44501 0.0116 | 

RA чы 


160 Degrees i | | ә NN 


















































47.1867 0.0073 | 0.4450 0.0033 
140 K Hz 35.3850 5.6726 | 0.0062 | 0.4450 0.0028 
| 44.5260 5.3757 0.0083 0.4450 0.0037 


29.2012 2.3823 | 0.0123 
1.5898 | 0.0126 
































16.2833 0.4450 0.0059 
| 
| 

3.1567 | 0.1545 | 0.0204 
Ши eee о 





170 Degrees 


39.1083 6.4486 | 0.0061 0.4450| | 0.0027 | 
13.6250 3.6166 | 0.0121 
18.3778 1.5898 | 0.0116 
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Appendix A: Experiment Data 100 Micron Sphere 


180 Degrees A V O 


Frequency | Scattered Incident | Ratio of Radius, Normalized 
Voltage Peak Voltage Peak Scattered & | | Voltage 


26.1592 2.3823 0.0110 
0.0112 
1.2265 0.0138 0.4450. 00061 - 
7.3292 0.8601 00065 


| 0.0174 1 0.4450 | 0.0078 
30 K Hz | 3.5700 | 0.1545 | 0.0231 | 0.4450 | 0.0103 



































190 Degrees 


















| 422913 6.4486 | 0.0066 10.4450 | 0.0029 
5.6726 | 0.0076 10.4450 0.0034 
$1122 0.0099 10.4450 
48.0690 | 3.6166 | 0.0133 0.4450 
90K Hz | 36.7250 |! 2.3823 | 0.0154 | 0.4450 0.0069 
80 K Hz | 22.3860 1.5898 | 0.0141 | 0.4450 
70KHz | 201035 1.2265 | 0.0164 | 0.4450 0.0073 
60 K Hz 10.3447 | 0.8601 | 0.0120 
50K Hz 9.9688 0.5325 | 0.0187 
40KHz | 6.4400 0.3127 0.0206 
30KHz | 3.6900 


| 






|S | 

200 Degrees | EN À 
150 K Hz 0.0073 0.4450 0.0032 
0.0102 0.4450 0.0045 
4.8513 | 0.0115 0.4450 0.0051 









100 K Hz | 41.5592 3.6166 | 0.0115 0.4450 0.0051 
90 K Hz 29.8485 2.3823 | 0.0125 0.4450 0.0056 






0.4350 0.0061 
0.0154 [04450] 0.0069 
0.0121 |04450| — 0.0054 
[ soKHz | 6.4110 | 055325 | 00120 [o 0.0054 
0.0177 
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Appendix A: Experiment Data 100 Micron Sphere 


210 Degrees | | 
Scattered | Incident | Ratio of Radius, Normalized 
Voltage Peak Voltage Peak | Scattered & Voltage 
15.0542 | 1.5898 0.0095 
70K Hz 16.3265 
50KHz 6.6985 0.5325 | 0.0126 


40K Hz | 5.0305 0.3127 | 0.0161 0.4450 0.0072 
30K Hz 0.0292 0.4450 0.0130 





| | NNNM 
101850 | 5. 0.0071 
417008 | 5. 2 0.0078 
46.0308 5. _ 0.0088 
100 K Hz 39.1815 | | 0.0108 
40 K Hz 
30K Hz | 2.9100 0.1545 
| | я 
[230 Degrees __ | 
130 K Hz 
120 K Hz 
100 K Hz 
90 K Hz 
80 K Hz 0.0147 


|. 0.0147 (0. 
70 K Hz 
60 K Hz 
_ 0.5325 
40K Hz | 5.8380 0.3127 0.0187 0.4450 
30K Hz | 44325 0.1545 0.0287 
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Appendix A: Experiment Data 100 Micron Sphere 











240 Degrees a AA A 
Frequency Scattered Incident Ratio of Radius, Normalized 












Voltage Peak | Voltage Peak Scattered & Voltage 
3.6166 — , _ 0.0115 (0.4450 
31.5850 2.3823 | 0.0133 
1.5898 | 0.0091 
1.2265 0.0136 0.4450 0.0061 
0.4140 0.4450 0.0071 
5.9080 | 0,5325 0.0111 0.4450 0.0049 
40 K Hz 3.6175 | 0.3127 | 0.0116 0.4450 0.0051 
30K Hz | 1.4740 | 0.1545 | 0.0095 0.4450 





кә 
UA 
© 
T 
e 
Ja 
- 
6 
e 
W 

















150KHz | 36.7973 ı 6.4486 | 0.0057 10.4450 | 0.0025 
140K Hz | 32.2718 | 5.6726 | 0.0057 | 0.4450 0.0025 
120 K Hz 43.9790 52102 | 0.0084 0.4450 0.0038 
110K Hz 41.8430 4.8513 0.0086 0.4450 0.0038 
100 K Hz 39.3550 3.6166 | 0.0109 0.4450 0.0048 
90 K Hz 29.6080 | 2.3823 | 0.0124 0.4450 0.0055 
80K Hz | 15.4360 1.5898 | 0.0097 0.4450 0.0043 
70 K Hz 11.0180 1.2265 | 0.0090 0.4450 0.0040 
60 K Hz 12.9610 | 0.8601 | 0.0151 0.4450 0.0067 

10.4450 | 0.0028 

260 Degrees | | 

150 KHz | 46.0317 6.4486 | 0.0071 0.4450 | 0.0032 








140 K Hz 48.6158 0.0086 0.4450| 0.0038 | 
11.3200 3.6166 0.0114 
14.6733 1.2265 0.0120 


55 










Appendix A: Experiment Data 100 Micron Sphere 





270 Пенев 


— о 
Voltage Peak — Peak а & Voltage 
36.7892 3.6166 0.0102 
60 K Hz 11.6650 0.8601 0.0136 0.4450 0.0060 
50K Hz 8.7938 | 0.5325 | 0.0165 
6.8925 | 0.3127 | 0.0220 0.4450 0. 0.0098 | 
0.0442 10.4450 | 
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Appendix B: Experiment Data 500 Micron Sphere 




















































Frequency Scattered Incident Ratio of | Normalized 

| Voltage Peak Voltage Peak Scattered & Voltage 

150 K Hz 95-3338 5.3727 0.0177 0.0080 

120 K Hz 72.9500 4.2008 0.0174 0.0078 

90 K Hz 34.8250 2.1966 0.0159 0.0071 
[ é0KHz | 163583 | 0387 | 0.0425  [04300| 0019] 
AAA | | —5 | 
е | | , | 





IS0K Hz | 84.5083 
29.5333 2.1966 0.0134 
0.3847 | 0.0456 0.4600 | 0.0210 











1 






110 Degrees 













IS0KHz | 94.4000 mn 4 0.0176 0.4530. 0.0080 

120 K Hz | 68.0333 4.2008 | 0.0162 10.4530 | 0.0073 

90KHz | 23.1333 2.4966 | 0.0105 | 0.4530 0.0048 

60 K Hz 10.8000 0.3847 | 0.0281 

30 K Hz 11.9300 | 0.0773 0.4530 
2o 2 









10.4550 | 0.0063 





ШШ |" 303500 : 2.1966 








60K Hz 15:296 7 0.3847 0.0398 0.4550 0.0181 





12.7250 0.1543 | 0.0825 0.4550 0.0375 


130 Degrees | 


0.4550 










0.0078 









0.4550 
NEM c _ 


140 Degrees 


0.0185 


























150K Hz | 72.3700 5:3227 | 0.0135 0.4570 0.0062 
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Аррепдіх В: 









Frequency 






150 K Hz 
120K Hz 
90K Hz 


160 Degrees 
150 K Hz 
120 K Hz 


60 K Hz 


30 K Hz 





120 K Hz 








180 Degrees 


150 K Hz 100.8000 5,3127 0.0188 0.4450 0.0083 


120 K Hz 


150 Degrees 





Scattered 
Voltage Peak 
62.8250 
84.8000 
36.4333 


85:3250 
83.9083 


90 K Hz 38.3750 


14.0833 
3.8133 


67.6000 


90.9667 





90K Hz | 


150 K Hz 
120 K Hz 
90 K Hz 









32.5583 


89.5250 
89.9917 
27.8833 











Experiment Data 500 Micron Sphere 


Incident 


3019927 
4.2008 
2.1966 


5.3727 
4.2008 
2.1966 
0.3847 
0.1543 








55972 


Voltage Peak 










| 















Ratio of 


Scattered & 


0.0117 
0.0202 
0.0166 


12.1000 0.3847 | | 0.00315 [0.4520 0.0142 
30K Hz | 4541 


0.0294 


0.0159 
0.0200 
0.0175 
0.0366 
0.0247 





90 K Hz 45.9917 2.1056” 0.0209 0.4550 0.0095 
60K Hz | 14.9250 0.3847 | 0.0388 | 0.4550 0.0177 
30 K Hz 3.3950 | 0.1543 | 0.0220 | 0.4550 0.0100 





60K Hz | 117917 0.3847 | 0.0307 0.4450 0.0136 
ЗОН 4.3500 0.1543 0.0282 0.4450 0.0125 


190 Degrees СИ 


0.0167 























Normalized 


Voltage 
0.4520 0.0053 
0.4520 0.0091 


0.0075 







0.4520 0.0133 


= 


! 0.4500 0.0111 

































0.0096 
04450| 0.0066 







0.4430 0.0074 


pU K Hz | 10.7333 0.3847 0.0279 0.4430 0.0124 
30 K Hz 3.2400 0.1543 0.0210 0.4430 0.0093 


— 
29 Degrees | | | |] 
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Appendix B: Experiment Data 500 Micron Sphere 






210 Degrees 
Frequency | 


























Scattered Incident Ratio of Radius, Normalized 
Voltage Peak Voltage Peak Scattered & Voltage 


11966 | 00137 10.4430] 00061 | 
















МЕНЕ |) | 
A ДЖ S SSS 
12008 | 00179 
33.0000 2.1966 | 0.0150 
60K Hz | 13.8337 .| 03847  : 0.0360 0.4450; 0.0160 





30K Hz 3.0500 0.1543 0.0198 | 0.4450 0.0088 











230 Degrees — — 
150 K Hz 88.4917 DSZ 0.0165 0.4450 0.0073 








120 K Hz 76.9917 4.2008 0.0183 0.4450 0.0082 
90 K Hz 42.0100 0.0191 0.4450 0.0085 


240 Degrees 
150 K Hz 89.5167 53727 | 0.0167 0.4490 0.0075 


BOX Hz 83.4667 | 4.2008 | 0.0199 0.4490 0.0089 
90 K Hz 36.0033 2.1966 | 0.0164 0.4490 0.0074 













60K Hz | 11.7692 | 0.3847 | 0.0306 0.4490 00113177 
30 K Hz 8.8000 | 0.1543 0.0570 ! 0.4490 0.0256 







ЕЕ. | - — 
R | 
105.2917 | 
79.0583 | ww 


90 K Hz 31.3250 2.1966 0.0143 0.4540 | 0.0065 


60 K Hz 11.9083 0.3847 0.0310 0.4540 0.0141 


30 K Hz 7.1800 | 0.1543 0.0465 0.4540 0.0211 


260 Degrees 


98.1417 5.3727 0.0183 
86.1750 4208 7 0.0205 
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Appendix B: Experiment Data 500 Micron Sphere 


270 Degrees 


Frequency Scattered Incident Ratio of m 
Voltage Peak Voltage Peak Scattered & Voltage 


150 K Hz 96.2083 5.3727 0.0179 0.4700 0.0084 
33.9167 2.1966 0.0154 
60 K Hz 14.3000 0.384 0.0372 0.4700 
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Аррепдіх С 
Input to Kargl's Program: 100 Micron Case 


Parameter | Name 10% Damping 


1000 1000 1000 

Bulk Modulus, K 2.25x10? 2.25x10? 2253600? 
Shear Modulus, N | em | оо | oo | 
kd 












162 (Y c4 





density, p 


porosity, P beta 0799999 0.99999 0.99999 


FATIH EF | чо | 
min wave number xlmin ЕЕС | 
max wave number | xımax | so | 


1.0 


















radius 0.332 
density, p Bulk 2231 
Bulk Modulus, K 3 5x10" 
Shear Modulus, 2.6x10°° 
porosity, P 0.31 
permeability, k 6. 5x19" 
à = (2 
1 ; D 


( 
1000 

Bulk Modul, K 2.25x10° 2.25x10° 
NEUE SS». | оо | 


[| = f] 
Bulk Modulus, K 2.25x10? 2.25x10? 2.25x10? 
0.001 


| dist | — os 2 | 
ee ЕНГЕН 1 
| nstert | o | 


1 

nend 75 
Lm | 1 
ENS | i 
4 


1.0 


ІШ 












`J 
(л 


Poe) | 
a 
s 
361 361 
Е |... 5 | PES 

30,000 _ 


OO fregmin| 30,000 30,000 30,000 
|. | déreg | 30,000 30,000 30, 000 
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Аррепдіх С 
Input Data for Kargl's Program: 500 Micron 


Parameter | Name 10% Damping 


1000 1000 
2.25x10° 225х107 
0,99999 0.99999 0.99999 


1.0 1.0 
max wave number 


— 0.332 0.332 
о 2231 2231 
Bulk Modulus, K 3.5x10' 3.5x10' 
Shear Modulus, 2.6x10°° 


STT 0.32 


574x107 
хе 0 g 
| density, p | mmo | 1000 | 1000 1000 
Bulk Modul, K | xe | зае | 2.2510 | ` 2.2520" 
viscosity, n | — 


1000 1000 1009 
Bulk Modulus, K 2.25x10? 2.25x10? 2.25x10? 
0.001 0.001 0.001 

















no Аа = 

density, f 
Bulk Modulus, K 
Shear Modulus, N 
porosity, P 
permeability, k 


min wave number хітіп 






361 361 





30,000 30,000 30,000 
30,000 


Аррепдіх С 
Input to Kargl's Progam: Rigid Solid Case 


Parameter Rigid Sphere 


0 
e 
= 
U 
) 


densit р 1000 


Bulk Modulus, K (-2.25x10?, 0.0)* 
Shear Modulus, N (-2.25x10?, 0.0)* 
porosity, P 0.999999 
permeability, k 

min wave number 

max wave number 


radius 0.332 
density, f 100,000 
Bulk Modulus, K (-2.25x10”, 0.0)* 
Shear Modulus, N (-8.44x10”, 0.0)* 
0.000001 


permeability, k као deci: 2 


(A) 
Q 
е 
ә 
е 
) 
) 


Це 


m 
e je - . 
о 


porosity, 


rhof 
kf 


density, p 1000 


Bulk Modul, K 2.25х10° 


viscosity, eta 


TO 100,000 
Bulk Modulus, K D 
viscosity, 0.01 


361 


nfreo 


`J 
іл 


*complex # format (real, imag) 


(л 


freqmin 30,000 


30,000 





dfreq 
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